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PREFACE 


This  report  presents  the  final  94  GHz  atmospheric  emission 
measurements  performed  during  February  and  March  of  1980  under  the  terms 
of  Contract  No-  N00173-780C-0165  between  the  Naval  Research  Laboratory 
(NRL)  and  the  Georgia  Institute  of  Technology  Engineering  Experiment 
Station  (GTEES).  Two  field  exercises  were  conducted  at  NRL  during  the 
program,  and  these  have  been  discussed  in  previous  reports.  The 
measurements  described  herein  were  obtained  with  a  slightly  modified 
version  of  the  radiometer  used  during  the  1978  and  1979  measurement 
programs.  Theoretical  models  of  atmospheric  fluctuations  are  briefly 
discussed.  The  data  taken  during  February -March,  1980  are  reviewed  and 
discussed,  with  chosen  data  analyzed  with  current  theoretical 
considerations.  Recommendations  and  conclusions  resulting  from  this 
work  are  presented. 


I.  INTRODUCTION 


The  turbulent  atmosphere  is  an  inhomogeneous  medium  in  which  the 
refractive  index  is  a  function  of  position  and  time.  Scintillation,  the 
observed  fluctuation  in  intensity  and  apparent  position  of  small  angular 
size,  is  caused  by  random  fluctuations  in  the  refractive  index  of  the 
earth's  atmosphere  through  which  the  signals  propagate.  The  majority  of 
the  investigations  of  atmospheric  turbulence  has  been  confined  to  the 
optical  and  the  low  frequency  regions.  Until  recently,  the  effects  of 
atmospheric  turbulence  on  millimeter  wave  propagation  have  received 
little  attention,  and,  in  addition,  most  millimeter  wave  investigations 
have  been  directed  to  one-way  link  transmissions  employing  coherent 
transmitters.  Only  a  few  studies  have  concentrated  on  passive 
radiometric  observations  of  atmospheric  fluctuation  effects. 

The  objective  of  these  upward  looking  atmospheric  measurements  was 
to  establish  the  spatial  and  temporal  characteristics  of  atmospheric 
absorption/emission  in  the  94  GHz  region.  This  atmospheric  information 
is  important  for  determining  limitations  on  radio  astronomy  resolving 
power,  and  applications  in  communications,  and  is,  in  turn,  necessary 
for  providing  the  overall  target  to  background  clutter  threshold  as 
viewed  from  space,  when  combined  with  high  resolution  data  on  target 
signatures  and  surface  clutter. 

In  this  report,  some  of  the  effects  expected  from  theory  are 
discussed  in  order  to  provide  a  preliminary  estimate  of  the  effects 
observed  in  the  experiments  performed  at  Georgia  Tech.  Most  of  the 
theoretical  coi  siderations  are  extensions  from  optical  studies;  this 
practice  is  questionable  until  further  data  and  analysis  exist  in  the 
millimeter  wavelength  region.  However,  some  millimeter  wave  radiometric 
studies  have  been  performed,  and  our  discussion  will  rely  heavily  on 
these  investigations. 
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Thp  effort';  of  at.rinsphoric  fluctuations  have  been  extensivoly 
sttifliof^  in  tho  nntirai  reninn  v/ith  little  consideration  qiven  to  the 
effects  in  the  nillineter  vfavelenqth  region.  Of  the  measurenents 
nerfome'f  at  microwave  frequencies  and  the  few  in  the  millimeter  region, 
practica''ly  a^l  have  been  for  propagat  ion  between  two  points  with 
coherent  sources.  bnly  a  fev/  have  treated  observations  or  relevant 
theory  for  nassive  raifiometric  observations.  '\s  a  result,  much  of  the 
theory  and  terminolooy  must  not  only  be  adanted  to  a  new  spectral  reoion 
blit  could  be  erroneously  used  to  internref  passive  observations.  I't  the 
current  stane  of  activities  in  the  area  of  interest  to  this  program, 
there  is  a  need  for  detailed  investinations  of  both  theory  and 
measurements.  For  the  few  investinations  of  atmospheric  fluctuations  by 
millimeter  wave  radiometry,  the  interest  has  been  in  slower  fluctuation 
rates  than  those  investigated  in  the  measurements  performed  in  this 
study. 


In  addition  to  the  above  considerations,  the  observed  fluctuations 
originate  from  different  causes.  In  the  visible  wavelength  region,  the 
refractive  index  variations  are  associated  with  the  fluctuation  of  the 
temperature  of  the  atmosphere  due  to  turbulence.  At  longer  radio 
wavelengths,  tte  variation  of  the  electron  density  in  the  ionosphere  is 
thp  cause  of  scintillation.  At  shorter  radio  wavelengths,  where  thp 
influence  of  the  ionosphere  is  diminishing,  variations  in  the  v/ater 
vapor  concentration  of  the  troposphere  become  increasingly  important  in 
causing  fluctuations  in  the  index  of  refraction.  Kemp  [11  has  pointed 
out  that,  in  all  of  these  regions,  the  medium  is  essentially  non- 
absorbing  and  attenuation  of  the  signal  in  the  turbulent  region  nlays  an 
insignificant  role.  For  millimeter  and  submillimeter  wavelengths, 
hov/pver,  this  is  no  longer  the  case.  Thus,  it  is  necessary  to  consider 
fluctuations  in  atmospheric  attenuation  in  order  to  determine  the  amount 
of  signal  fluctuation  present  in  a  '"adiometric  observation.  In  most  of 
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the  analysis  performed  in  this  spectral  region  [2J,  the  model  uses  tlie 
properties  of  the  dry  atmosphere  to  conform  to  a  standard  atmosphere, 
e.g.,  the  US  Standard  Atmosphere  [3],  The  water  vapor  concentration  is 
assumed  to  have  a  mean  profile  but  the  concentration  is  allowed  to 
fluctuate  from  point  to  point  around  this  mean  profile. 

One  of  the  earliest  investigations  of  passive  observations  of  the 
fluctuation  components  of  the  atmospheric  noise  temperature  was 
performed  by  Orhaug  [4]  at  P.  GHz.  He  employed  a  12’  parabolic  antennn 
at  NRAO  and  considered  the  small  scale  fluctuations  in  brightness 
temperature  during  periods  with  no  precipitation  to  be  an  important 
limitation  in  radiometry.  His  short  fluctuations  were  on  the  order  of 
1-2  ”k  with  an  average  periodicity  on  the  order  of  a  few  minutes.  The 
system  integration  time  was  5  seconds,  long  compared  to  our 
requirements.  In  Orhaug's  observations  [4],  the  atmospheric  effects 
which  he  described  were  due  to  the  thennal  emission  by  the  atmosphere, 
resulting  from  absorption  characteristics  of  the  transmission  medium. 
The  influence  on  phase  characteristics  of  a  wave  propagating  through  the 
medium  was  not  included.  Orhaug  has  used  results  published  by  Hogg  [5] 
to  demonstrate  the  effects  of  changing  water  vapor  content  within  the 
receiver  beam  as  a  function  of  zenith  angle.  Hogg  used  the  brightness 
noise  temperature  from  a  standard  atmosphere  having  a  water  vapor 
content  of  10  g/m  at  ground  level  and  from  a  humid  atmosphere  with  20 

3 

g/m  .  In  Figure  1,  the  effect  of  increasing  the  water  vapor  content 

3 

from  10  to  20  g/m  is  indicated  for  10  GHz.  At  the  zenith  direction, 
it  is  seen  that  the  antenna  noise  temperature  difference  for  the  two 
water-vapor  contents  is  4°K  increasing  with  zenith  angle.  The  effect 
of  different  water  vapor  contents  in  the  propagation  path  increases 
considerably  as  frequency  is  increased  into  the  millimeter  region. 

A  very  rough,  but  interesting  estimate  of  the  variation  in 
brightness  temperature  can  be  made  for  the  assumption  that  a  change  in 
absorption  coefficient  occurs  over  a  limited  height  interval  aL. 

at  =  Te'“^^(AaAL)  (1) 
where  the  absorption  coefficient  is  «  +  Au  in  the  interval  aL. 
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In  this  case,  Aa  is  a  function  of  pressure,  temperature  and 
water-vapor  variation,  which  is  the  most  significant  parameter.  Orhaug 
has  considered  the  results  of  refractometer  measurements  to  obtain  the 
variation  of  the  index  of  refraction  as  a  function  of  the  variation  of 
the  water  vapor  content:  The  index  of  refraction,  n,  is 

„  .  ,  .  ,Z5p  .me  ,  ,36 

T  T  j2 

where  P  is  the  total  pressure  in  mb  and  e  is  the  partial  water-vapor 
pressure.  With  the  refracti vity,  N  =  (n-1)  10^,  it  is  assumed  that 
variation  in  N  is  caused  only  by  variation  in  e,  then 


an  =  4.2Ae  for  T  =  BOO^K 


Orhaug  indicates  that  refractometer  measurements  at  different  altitudes 
have  given  large-scale  variations  in  N  of  the  order  of  50-200  N-units 
for  the  first  5  km  of  the  earth's  atmosphere.  For  AN=200,  the 
corresponding  e-variation  is  ao=50  mb. 

The  absorption  coefficient  for  water  vapor  is  directly  proportional 

3 

to  the  water  vapor  density,  0  grams/  m  and 


So  the  corresponding  change  in  brightness  temperature  is 

aT^  Te'^n^''4a^^-AL)  (4} 

where  is  the  absorption  coefficient  in  nepers.  Orhaug  has  an  error 

of  a  factor  of  10  in  his  calculation  so  that  the  fluctuation  AT.  is 

b 

not  as  large  as  he  indicates.  In  our  case,  with  the  correction  made, 
aT^  =  1  .75°K 

fora=0.3  dB/km  =  .0346  neper/km 
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.1  =  5  gm/m^ 

Ae  "  r.  o 

T  =  250  K  e  =  0-2 

0  =  600  rnn  Hg  aL  -  1km 

This  variation  in  the  brightness  temperature  is  observable  with  our 
system,  but  is  not  as  large  as  was  observed  in  some  measurements. 

Sollncr  [Gj  has  made  measurements  of  the  frequency  spectrum  of 
fluctuations  in  submillimeter  sky  emission  and  absorption.  The 
observations  were  for  single  beam  techniques  and  double  beam  techniques 
(at  fixed  separation  at  6  feet).  The  results  were  analyzed  in  terms  of 

_3 

their  power  spectral  density  between  4  x  10  Hz  and  0.25  Hz.  Tor 
observations  involving  the  fluctuating  atmosphere,  Sollner  has 
considered  the  observed  brightness  of  a  source  of  brightness  B^(  v  ) 
propagating  through  an  emitting  and  absorbing  medium  of  brightness  B^(v) 
and  optical  depth  t(v).  This  observed  value  is 

B  =  yB^(v)f(v)e'^^'')dv  +  (v)f(v)  (1-e-"("))dv  (5) 

where  f(v)  is  the  normalized  response  of  the  detector;  the  first  tenn  is 
the  contribution  from  source  (transmission  term);  the  second  term 
results  from  the  intervening  medium  (emission  term).  The  fluctuations 
in  these  terms  are  transmission  noise  and  emission  noise.  If  one 
assumes  an  effective  optical  depth,  the  transmission  term  can  be  written 

B^  =  e'^ej  B^(v)f(v)dv  =  e'^e  (6) 

Tor  a  source  constant  in  time,  any  fluctuations  in  the  transmission 
term  are  due  to  change  in  This  term  can  be  made  to  dominate  the 

omission  noise  by  choosing  a  sufficiently  bright  source.  Sollner  gives 
T  as  +  ATg(t)  where  Atg(t)  is  the  fluctuations  of  the  effective 


f 

I 

I 
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optical  depth 


Bj  =  B^exp  -[  Tg  +  '.ig  ]  (7) 

and  the  power  'ipectrum  of  as 

S.  _(f)  =  It,  f"'  '«) 

'e 

where  the  riapnitudo  of  depends  oni^f  is  the  frequency  in  Hz  and  i,  is 
a  constant  detemined  from  observations. 

In  the  case  of  emission  noise,  simplifications  are  not  possible 
because  changes  v/ith  time  as  does  ^ .  The  physical  situation  of 
regions  of  differing  and t  passing  through  the  observed  solid  angle 
can  be  represented  in  terms  of  spatial  and  temporal  distributions  for 
and  i.  In  the  investigations  performed  in  this  program,  neither  a 
tracking  capability  nor  sufficient  resolution  to  localize  solar  areas  of 
constant  brightness  were  available  so  that  the  important  aspects  of  the 
absorption  term  could  not  be  investigated.  Such  observations  should  bo 
performed  in  the  future.  Further  work  should  be  performed  in  other 
well-defined  observing  windows  and  correlations  should  be  made  vnth 
temperature,  water-vapor  and  wind  velocity. 

Several  additional  considerations  can  be  given  to  millimeter 
fluctuation  effects.  In  the  case  of  millimeter  wave  propagation 
observations,  theory  predicts  a  strong  dependence  of  the  scintillation 
amplitude  and  angle  of  arrival  variations  on  the  himidity  structure 
parameter  ^  in  addition  to  the  temperature  structure  parameter  Cj. 

Tatarski  [7]  has  shown  that  the  resulting  energy  distribution  in 

the  turbulent  atmosphere  is  log  normal,  characterized  by  a  variance 
2 

ic  that  is  a  function  of  the  degree  of  atmospheric  turbulence. 
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Chornov  [R]  and  Tatarski  [7],  in  their  original  work,  treated 
mainly  optical  fluctuations  and  neglected  the  effects  of  absorption  on 
the  fluctuations.  Recent  work  of  Russian  workers  has  considered 
fluctuations  in  the  millimeter  and  submi 1 1 ineter  wavelength  regions, 
requiring  the  inclusion  of  absorption  by  atmospheric  water  vapor. 
l7yumov  rg]  has  solved  the  wave  equation  to  account  for  absorption 
resultinn  in  expressions  for  amplitude  and  phase  fluctuations  valid  for 
millimeter  wave  propagation.  As  a  result  of  this  t«rk,  the  index  of 
refraction  N  is  given  by 

M  -  n  +  im 

with  n  =  n  ^  u  19) 

o 

and  m  =  m  +  v 
0 

Mere,  n^^  and  are  mean  values  of  the  real  and  imaginary  parts  of  N, 

and  and  V  are  the  fluctuating  parts. 


Armand  [1(1]  has  given  the  spectra  of  fluctuations  of  the  real  and 
imaginary  parts  of  the  index  of  refraction  and  their  cross-correlation 
in  terms  of  the  temperature  and  humidity  fluctuations  [See  McMillan  et 
al ,  Reference  11]. 


Hurvich  [12]  has  given  values  of u  and  vfor calculations  of  the 
spectra  of  fluctuations 


;.  =  (  Ki^  +  K2  y  +  )  x  1 


-6 


r(  P  ,  T  ,  e  )  ^  -1^  y-  X  10'^ 

0  0  0  e  T  2  P„ 

0  '  0 


(10) 

(11) 


=  yS^K/mb 

=  72'’K/mb 

=  3.7  X  10^  (°K)^  (only  weakly  dependent  on  \) 
p  =  atmospheric  pressure  in  mb 


e  =  partial  pressure  of  water  vapor  in  mb 

A  =  transmitted  radiation  wavelength 
e^,  p^,  Tg  =  stationary  values  of  e,  p,  T 
Y  =  absorption  coefficient  in  nepers/km 
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The  forms  of  the  spectral  distributions  of  fluctuations  of 
temperature  Tj(q)  and  humidity  fp(q)  are  qiven  by  burvich  [1?]  (see 
[11]).  For  the  relationship  between  C-j ,  the  temperature  structure 
parameter,  and  the  index  of  refraction  structure  parameter, 

tIcMillan  et  al  [11]  have  used 

=  C^(T^/79p)  X  10'^  (12) 

Actually  C  ^  and  are  the  meainqful  parameters,  but  several 

authors  quote  values  for  and  Cy. 

The  limitations  imposed  by  atmospheric  fluctuations  on  the  maximum 
linear  dimensions  of  larqe  telescopes  have  been  considered  by  Bastin 
[13].  The  limitation  of  angular  resolution  arises  from  differential 
phase  change  in  radiation  reaching  either  side  of  a  large  telescope  as  a 
result  of  changes  in  refractive  index  of  air  along  the  extreme  rays. 
Bastin  has  considered  that  relatively  small  changes  in  total  water-vapor 
in  the  solar  direction  can  be  determined  as  a  function  of  time  from 
small  fractional  variations  in  the  transmitted  solar  intensity.  This  is 
the  consideration  that  Sollner  had  made  [6].  For  deducing  the 
fluctuation  effects,  it  is  assumed  that  a  fixed  distribution  exists  with 
random  fluctuations  in  the  concentration  of  water-vapor  with  respect  to 
the  atmosphere  drifts  through  the  antenna  beam  due  to  wind  movement. 
This  widely  used  assumption  is  in  some  cases  referred  to  as  the  Taylor 
hypothesis,  but  is  actually  Tatarski's  hypothesis  of  "Frozen-in" 
turbulence.  Under  this  assumption,  temporal  changes  can  be  related  to 
spatial  ones  and  therefore  to  phase  changes  which  would  be  expected  to 
occur  between  the  spatial  limits  of  a  large  telescope.  Kemp  [2]  has 
extended  this  work  by  employing  concepts  put  forth  by  Brooker  [14]. 

Kemp  [2]  has  given  the  absorption  and  refractive  index  of  water- 
vapor  asGi(v)  and  n(\))  by  the  relations- 

'i{\>)  =  K(v)p  =  absorption  coefficient 
n(\>)  -  1  -  L(v)p,  n(v)  =  refractive  index 
and  p  =  density  of  water  vapor 
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Ho  assimes  that  all  i rrequl arities  have  (''^sontially  the  sane 
Himension  D  correspondi nq  to  the  scale  of  larqest  significant  eddies  in 
the  turbulence  and  uses  for  the  variation  in  water-vapor  density  the 
mean  square  value  for  the  deviation  from  the  mean  ('f))'  .  This 
variation  produces  a  corresponding  change  in  refractive  index  and 
attenuation  of  the  signal  within  the  region. 


The  fluctuation  in  refractive  index  is  given  by 

The  change  in  signal  intensity  is 

.",1  =  I  A I  for  At  ■  '  1 

J  .1  V  V 

where  L.  =  signal  intensity 
At  =  K(v)  DAo 


(13) 


Thus, 


A I 


change  in  optical  depth  of  the  region 

K(v)n  (ap)^ 


(14) 


=  relative  mean  square  intensity  fluctuation  per  region 

The  mean  square  fluctuation  in  phase  of  a  wave  of  wavelength  A 
propagating  through  one  region  resulting  from  a  refractive  index  change 


IS 


(A^,)|  =  (2^)  (An)^ 


2  n  L(  v) 

rr  RT^. 


-12  TAI  l2 

J  M 


AI 


(15) 


where  S(j) 


2tivL(  v) 

cKTvy 


=  the  scintillation  coefficient 
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[f  now  a  wave  travels  a  distance  through  the  fluctuation  layer,  it 
will  encounter  regions.  With  the  assimption  that  the  v/ater 

fluctuations  are  uncorrelated,  the  phase  fluctuations  will  add  in  random 
walk  fashion 


(16) 


Similarly,  the  relative  fluctuations  in  intensity  are 


(■A 


1  /  D 

\' 


2l  /  V' 

'v/' 


(17) 


Thus,  the  phase  fluctuation  for  a  v/avc  passing  through  a  layer  is 
related  to  the  intensity  fluctuation  by 


(^!)^  =  S^(v) 


(18) 


for  (  1  .  The  scintillation  coefficient  S(v)  is  inversely  proportional 
to  pressure  and  is  a  weak  function  of  temperature  and  the  concentration 
of  water  vapor.  The  coefficient  can  be  taken  as  a  constant  throughout 
the  troposphere,  and  the  mean  square  fluctuation  in  signal  phase  can  be 
determined  from  intensity  fluctuation. 


Kemp  [2]  has  considered  the  practice  of  determining  the 
fluctuation  in  phase  difference  for  the  wave  arriving  at  two  points 
separated  by  a  distance  d,  which  is  perpendicular  to  the  wave  being 
propagated.  If  C(d)  is  the  correlation  coefficient  between  the  phase 
fluctuations  at  the  two  points,  the  phase  difference  between  the  signals 
at  the  two  points  is 

-  4>2)1^  =  2{l-C(d)}  (A*)^  (19) 

In  addition,  let  =  the  angular  scintillation,  i.e.  the 
fluctuation  in  the  direction  of  arrival  of  the  phase  front. 


T 


(20) 
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To  consider  the  atmospheric  limits  on  an  instnnent  of  aperture  d, 
Kemp  compared  i^with  the  minimum  resolved  angle  due  to  diffraction.  For 
the  instrument  to  be  atmospherically  limited, 

"s  "d 

or  2  jl-C(d)|  (A-lO^  •  (2.44.)^  (21) 

Thus,  when  the  random  phase  fluctuation  due  to  the  medium  exceeds 
1.74-'  radians,  for  some  value  of  d,  the  atmospheric  angular 
scintillation  will  exceed  the  diffraction  limit  so  that  the  instrument 
is  atmosphere  limited. 

2 

If  the  mean  square  difference  coefficient  A  (d)  is  related  to  the 
correlation  coefficient  by 


A^d)  =  2jl-C(d)J  , 

(22) 

then, 

A^(d)  •  (Ao)^  -  (2.44,,)^ 

(23) 

In  the  millimeter  wavelength  region,  phase  fluctuations  and  the 
difference  coefficients  are  not  measured  directly,  but  they  can  be 
determined  from  observable  intensity  fluctuations  due  to  the  variation 
in  atmospheric  attenuation. 

For  1^-  1,  the  variation  in  phase  and  intensity  are  directly  related 
to  fluctuations  in  water-vapor  density.  Therefore,  the  correlation 
coefficient  and  the  mean  square  difference  coefficient  for  phase  are 
identical  to  those  for  intensity  fluctuations. 
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With  equations  IB  and  23 


A I 

'."(d)  S^(v)  1-j-l'  (2.44n)' 


where  now  (d)  is  for  intensity  fluctuations. 


The  instrument  is  atmospherically  limited  if 


(24) 


.'.1 

V 


(1.74  )^/S^(v). 


(25) 


In  addition,  the  limiting  size  of  a  diffraction  limited  instrument 
can  be  obtained  from 


A^(d)  =  (1  .74Tr)^  S’^(v)  j 


A I 


I 


•1 


(26) 


Neither  the  work  of  Kemp  nor  that  performed  in  this  program  could 
measure  intensity  fluctuations  from  two  positions  but  it  has  been 
possible  to  measure  intensity  fluctuation  as  a  function  of  time  at  a 
given  position.  It  is  therefore  possible  to  use  Taylor's  hypothesis  to 
relate  the  mean  square  difference  coefficient  for  a  time  interval  t  to 
the  mean  square  difference  coefficient  for  the  spatial  separation,  d,  by 
the  relation 


2  2 

=  o(t)  =  A^  =  o(d) 


(27) 


t  =  V  ’  ”  drift  velocity  of  the  atmosphere 

Kemp's  v/ork  indicated  that,  during  his  measurements,  the  drift  time 
interval  for  fluctuations  was  about  40  seconds,  which  when  combined  with 
the  wind  speed  of  10  m/s  gives  the  turbulent  region  dimension  to  be 
400  m. 
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In  their  i nvestiqations  of  millimeter  wave  atmospheric  fluctuation 
effects,  several  authors  have  employed  various  concepts  and  methods  of 
expressinc]  the  effects.  As  previously  indicated,  the  majority  f)f 
treatments  applies  to  active  one-vwy  propagation  and  many  discussions 
are  not  applicable  to  passive  observations.  It  is  important  for  this 
discussion,  however,  to  include  brief  comments  on  these  various 
i nvestigations: 


1.)  The  Kolmogorov  model  [15]  assumes  homogeneous  and  isotropic 
conditions  of  the  atmosphere  to  describe  the  index  variations.  For  a 
particular  range  of  separation  between  two  points,  and  r^*,  the 

model  yields 


n{r^)  -n(r2) 


(28) 


where  denotes  an  ensemble  average  and  is  the  index  structure 
constant.  The  separation  range  for  validity  of  the  model,  often  referred 
to  as  the  inertial  subrange,  is 


? 

0 


L 


0 


(29) 


v/here  L  and  are  the  outer  and  inner  scales  of  turbulence 

respectively.  and  may  be  thought  of  as  the  approximate  maximum  and 

minimum  of  the  eddy  size.  In  the  atmosphere, ranges  from  a  millimeter 

to  centimeters,  whereas  for  horizontal  propagation  in  the  low 

atmosphere,  is  about  1/3  the  height  above  ground.  For  separations 

greater  than  L  ,  the  mean  square  index  fluctuation  levels  off  to 
2  2/3  ° 

C  v/hereas,  for  separations  less  than  viscosity  effects  cause  a 

very  rapid  decrease  in  index  fluctuations. 


14 


2. )  For  intensity  fluctuations  in  the  millimeter  wave  region, 
theory  requires  consideration  of  the  problem  in  two  separate  domains 
dependent  upon  the  si/e  of  the  outer  scale  of  turbulence  compared  to  the 
first  Fresnel  zone  along  the  propagation  path  of  length  K.  The  cases 
are 

L  '  AR  and  L  ■  /aR" 

0  0 

Most  rough  estimates  of  turbulence  effects  in  the  mm  wavelength  region 
are  based  on  Tatarski's  calculations  [7J  ,  valid  for  ■  Ar  .  flo 
simple  quantitative  models  of  amplitude  fluctuations  for  the  case  of 

L  '• 

0 

exist,  but  Tatarski  for  a  plane  wave  with  L  >  A¥  gives  the  variance 

0 

of  the  log-intensity  fluctuations  as 

.2  =  <  (10  log^Q  i  > 

=  23.39  (30) 

where  k  =  2i/> 

Worst-case  estimates  [16]  for  C  have  been  made  on  basis  of 

2  ^13  -2/3 

optical  measurements  to  give  C  6  x  10  m  for  strong  turbulence. 

3. )  Using  optical  constants  for  obtaining  estimates  of  millimeter 

wave  turbulence  effects,  one  must  realize  that  optical  turbulence  is 

mainly  dependent  upon  atmospheric  temperature  fluctuations  and  that 

varying  water  vapor  effects  are  negligible.  For  millimeter  waves,  water 

vapor  contributions  to  the  index  of  refraction  become  important.  Brown 

has  shown  that,  for  microwaves  (>  10  GHz)  statistical  variations  in 

2 

water  vapor  below  8  km  can  produce  values  of  C  that  are  more  than 

two  orders  of  magnitude  greater  than  the  values  for  the  corresponding 

2 

optical  case;  therefore,  estimates  ofo  based  on  optical  constants  may 
be  in  serious  error. 
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4.)  Armand  et  al  [16]  examined  fluctuation  effects  near  the  920 
.:'Ti  water  line  and  found  that  on  the  absorption  line  center,  the 
amplitude  fluctuations  were  approximately  five  times  less  than  in  the 
980  II m  wi  ndow. 


6.)  Mavrokoukoul aki s  et  al  [17]  have  compared  the  measured 

variances  of  log  amplitude  fluctuations  at  3G  and  110  GHz  as  a  function 

of  time  and  showed  that  the  fluctuations  at  the  different  frequencies 

were  very  well  correlated.  Ho  et  al  [18]  have  performed  simultaneous  mm 

2 

and  X-band  refractivity  measurements  of  C  obtaining  over  a  one-hour 

2  "  -1  a  _?/2 

period,  respective  average  values  of  C  of  0.25x10  m  and  0.32 

^  ,n-14  -2/3  " 

X  10  m  . 

6.)  Andreyev  et  al  [19]  have  made  measurements  at  A  =  2  mm  on  a 
horizontal  path  of  length  5.6  km  in  a  strongly  turbulent  atmosphere. 
They  were  concerned  about  testing  Tatarski's  hypothesis  of  "frozen-in" 
turbulence  and  estimated  the  width  of  the  spectrum  on  intensity 
fluctuations  due  to  cross-transfer  of  homogeneities  across  the  path  of 
propagation  to  be 

'  F  =  <-  v^>//^'  131) 

where  <  Vj^  =  the  mean  speed  of  homogeneities  transfer 
and  L  =  path  1  ength 

Andreyev  et  al  found  aF  =  0.17  Hz.  They  divided  the  fluctuations 
into  fast  (■  0.1  Hz)  and  slow  (  <  0.1  Hz).  Their  conclusion  v/as  that 
"frozen-in"  turbulence  doesn't  describe  intensity  fluctuation  quite 
correctly.  The  conclusion  was  that  "frozen-in"  turbulence  was  confirmed 
when  it  describes  intensity  fast-fluctuations.  Intensity  slow- 
fluctuations  are  assumed  to  be  caused  by  cross-transfer  and  evolution  of 
large  scale  inhomogeneities  whose  sizes  are  more  than  AU  .  To  estimate 
Cp,  Andreyev  et  al  considered  the  effect  of  aperture  averaging  and 
obtained  values  for  C^in  the  range  from  0.12x10"^  to  0.61xl0‘^m"^^^. 
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In  support  of  measurement  of  atmospheric  turbulence,  a  vertical 

profile  of  the  thermal  structure  of  the  atmosphere  would  be  very 

important.  Button  [20]  combined  thermal  sensor  technology  for 

microthermal  measurements  with  radiosonde  balloon  systems.  This 

resulted  in  an  extension  of  turbulence  sensing  to  heights  up  to  2b  km 

2 

above  sea  level.  This  measurement  technique  provides  Cy  (h)  data  where 
2 

Cy  (h)  is  the  temperature  structure  coefficient  and  h  is  the  altitude. 

2 

The  refractive  index  structure  coefficient  C  (h)  is  obtained  from 

2  ^ 

relationships  with  (h)  at  least  for  optical  effects.  Bufton  obtained 
the  mean-square  temperature  difference  between  two  inicrothennal  probes 
as  a  function  of  altitude.  This  is,  by  definition,  the  temperature 


structure  function,  D^,  at  probe  locations  r^  and  r^: 

'^2^  "  -  T(r2)]^> 


T(r^)  =  temperature  at  point  r^ 

r  =  kj  -  r^l 


C-j-  is  a  strength  parameter.  A  larger  value  indicates  more 
temperature  fluctuations,  which  are  associated  with  more-turbulent 
mixing  of  air.  The  general  expression  for  refractive  index  as  a  function 

O 

of  temperature  and  wavelength  provides  the  connection  between  C-.  and 

2  ' 

_ _ _ 4-i__  _  r\ 


.  Bufton  uses  the  expression  2 

c2(h)  =  [79-9  P(h)  X  lO'^l  C^(h)  (32] 

where  P(h)  =  atmospheric  pressure  (mbars) 

T(h)  =  ambient  temperature  (°K) 
h  =  altitude. 

This  relation  has  been  used  by  McMillan  et  al  [11]  [See  Equation 
12],  but  Bufton  has  indicated  that  the  relation  applies  for  0.5  ym. 
Therefore,  caution  should  be  exercised  in  using  it  at  millimeter 
wavelengths. 


1) 


K'  - 


/'  rpcpnt.  piihl  i  rati  on  hy  Hill,  Clifford  and  Lawrpnco  [??]  trpats  the 
pf^'prts  of  re^'rartivp-index  and  ahsorption  fluctuation?;  on  propanation 
for  the  nirrowave  rpqion  throuoh  the  IR.  This  work  ran  ho  adapted  to 
thp  vertical  ra'^iomptric  ohservat i ons  of  interest  to  this  pronran.  The 
r'o,".''  nroup  das  i  nvpst  i  rjatpH  the  Hependenre  of  fluctuations  in 
atinspheric  ahsorption  and  refraction  upon  fluctuations  in  temperature, 
hiiniditv  and  pressure. 

Tor  the  applications  of  interest  in  this  report,  fluctuations  in 
atmospheric  refraction  are  not  significant.  The  work  of  dill  et  al  [P?! 
^as  considered  the  contrihutions  from  line  ahsorption  hy  H^n,  They 

L 

have  'ipvploppd  functions,  reiatim  the  fluctuations,  v^hich  arc  necessary 
for  evaluatinn  deqradation  of  el ectronaqnetic  radiationhy  turbulence. 
In  the  calculation  of  the  turbulence  effect,  v/e  must  choose  a  set  of 
mean  atnnspheric  conditions.  Tince  the  observations  are  for  vertical 
absorption  effects,  the  mean  atmospheric  conditions  are  chosen  for 
vertical  layers  of  the  atmosphere. 

In  reference  fP?],  it  is  assumed  that  turbulent  fluctuations  in 
total  pressure  oive  a  neqlioible  contribution  to  absorption  and 
refraction  fluctuations.  This  assumption  is  for  horizontal  active 
system  nronaqatinn,  but  is  assumed  to  apply  to  vertical  observations, 
bbether  this  is  a  reasonable  assumption  for  vertical  observations 
remains  to  bp  determined. 

In  the  millimeter  wavelenqth  reqion,  hiriidity  fluctuations  dominate 
absorption  nat.ions.  In  order  to  determine  the  effects  of 

absorption  fluctuations  on  antenna  temperature  IaT),  it  is  possible  to 
emnloy  the  variation  of  atmospheric  absorption  as  a  function  of 
fluctuations  in  temperature,  humidity  and  pressure.  It  is  necessary  to 
provide  these  functions  that  relate  the  fluctuations  in  order  to 
evaluate  the  deqradationof  radiation  by  turbulence.  It  is  necessary  to 
choose  a  set  of  mean  atmospheric  conditions. 
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Total  absorption  is  the  sum  over  all  lines- 
^  L.  P. 

n  j  -  i;/4ii V 

Q  =  absolute  humidity 
The  fluctuations  in  are  given  by- 


T  oS.  E.'-  V.  exp(-v./CT) 

ST  TT  "  FT  ■  CT  1  -  exp(-v./CT) 

I  6(1  . 

- =  -b  +  R 

!J  .  I 

1 

P  6a  . 

;:tV  =  1  « 

Q  6a  . 
a^  6(3 

where  R  =  — — 

P  +  4kQT 
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"i  -  j 

.  ...  g 

and  A  =  212  ,  b  =  0.62 


If  we  were  only  concerned  with  fluctuations  due  to  humidity 
fluctuations,  the  result  would  be 


> 


=  1  + 


=  1  + 


11 

<Q> 

'‘-i 

Q 

(3.  6u. 

a.  5Q 

a.  6g 

4kOT 

T  ^ 

P  +  4kQT 

<Q> 
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The  total  fluctuation  in  atmospheric  absorption  is 


-s  =  Z  “i 


6  q_ 

<Q> 


1  + 


\  g  V 


4kQT 
,P  +  4kQT 


) 


15 

<q> 
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In  the  low  frequency  Unit,  ■  CT ,  applyinq  to  the 

ni  1 1  irneter/submi  1 1  imeter  wavelength  region. 


S.(T)  exp(  -e|  /  CT) 

The  dependence  of  linewidth  on  P,  D,  T  is  given  as  [22] 


with  b  =  0.62 


=  width  of  line  i  at  reference  pressure  P  and 
temperature  T^.  ” 


Hill  et  dl  used  formulas  for  differential  changes  to  find  the 
fluctuations  caused  by  turbulence  (accurate  to  first  order  in  fl  uctuations) . 
Variables  P,  T  and  0  are  written  as  the  sun  of  their  mean  values  •'  P  '  , 

!•  and  Q  and  their  fluctuations  caused  by  turbulence  - 


P  =  -P:-  +  SP 
T  =  •'T.  +  ST 
Q  =  -Q>  +  SQ 


The  following  relations  [22]  are  applicable  to  the  measureinents  of 
this  investigation  - 

Imaginary  part  of  refractive  index  attributable  to  a  single 

absorption  line  =  n., 

1 1 

Absorption  coefficient  due  to  the  same  line  =  B.j 

n^j  = 


g  =  line-shape  factor 
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'ahio  1  of  l<pf  oronco  r?rl  nivo*.  tho  lino  shapp*",  anH  thpir 
•'tn-ivat  i  VPS  with  rpsppct  to  linowi-tth.  It  is  nossiblf'  to  pripl  oy  thp 
a'lovo  formulas  to  dptominn  thp  tl  urtuat  i  on  in  thp  absorption 
r nn f i r i rn t  nvnr  anv  vprtical  path  or  spofipnf  of  thp  path  as  a  funrtion 
of  thp  1 1  UP  t  iiat  i  ons  of  thp  threp  variahlos  ( F’ ,  T,  0).  Two 
annroxinat i ons  of  thp  fluctuation  in  hriohtnoss  tnnppraturp  havo 
rpppnti''  hopn  nnploypd  hy  P.  I!.  '^cMillan  of  our  laboratory  to  cornparp 
v/ith  thp  data  rpportpd  hnrn.  bp  ronsidprpd  fluctuations  in  thp  vprtiral 
.(atpr  vapor  distribution  to  nrodiicp  fluctuations  in  thp  absorption 
ropb^^iriont  watpr  vanor.  In  onp  case,  hp  considorod  a  fluctuation  in 
huniditv  at  oroiind  IpvpI  and  a  correspondi no  fluctuation  vertically 
tbroiinh  the  atnosphprp.  '’’his  approach  nave  fluctuations  in  brightness 
•■pnnpraturp  on  the  order  of  but  loss  than  the  reasiirpd  values.  Thp  use 
o'"  a  fluctuation  throuohout  the  atmosphere  is  not  the  most  reasonable 
assiinpt  ion.  ''  nore  reasonable  one  is  to  assijn’e  that  the  turbulence 
ornirrod  in  a  layer  approximate! v  1  '<n  vnde  at  some  altitude  fused  as  10 
kn  in  the  case  analyredl.  The  iustification  of  a  stratification  like 
this  night  bp  found  in  thp  work  of  Oufton  [I’ll.  The  calculations 
msiilted  in  brinhtness  tenperaturr  fluctuations  on  the  order  of  the  RtlS 
values  of  fluctuation  observed  in  the  experiments,  but  almost  an  order 
or  naonitijde  smaller  than  the  peak-to-peak  fluctuations  vfhich  have  been 
observed.  Pos'^iblv  a  stratification  of  the  turbulance  in  the  fom  of 
lavprs  or  horizontal  sheets  of  different  amplitudes,  compatible  with 
thermosondp  turbulence  data  would  be  the  most  applicable 

assumption. 
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II I-A. 


I’AniOMFTRir  NTASlIRf-MFIIT  ‘^YSKM 


Th.'>  dti^n'^^nhpric  'rpasiirennnt s  conductpd  durinn  this  contract  have 
tcpn  nhtainp'i  with  a  sunprhpterodvnp ,  double  sideband,  nicko-swi  f.chpd , 
cti7  radionptcr.  The  radioneter,  prpviotisly  dpscrihed  in  the  Serii- 
Annual  Ppport  is  dppictpd  in  f  iqiirp  ?.  Thp  nost  unusual  fpaturr 

of  thr>  radiopptpr  front-ond  i^  thp  puasi-optical  antenna  feed  which 
iilnv'/s  one  of  several  antenna  schemes  to  be  emploved.  Fioure  ?  shows 
the  front-end  as  it  was  oneratpd  at  the  prime  focus  of  the  Naval 
'•psearrh  Laboratorv's  IP  ^oot  dish  durinq  the  ]07P  and  1^70 

npasurp'’'pnt  nroqrams.  for  the  l^PP  neasurenents ,  the  radiometer  antenna 
feed  v/as  intenratpd  with  a  '’A  inch  Cassenrain  antenna  system  by 
rnnnertinq  the  conical  horn's  wavequide  port  to  thp  feed  horn  of  the 
casspqrain  antenna.  Antenna  patterns  for  the  modified  system  are  shown 
i n  Fi qures  and  A  , 

A  new  sharpless  wafer  mixer  vws  installed  in  the  radiometer  for  the 

loPO  measurements  v/hich  greatly  improved  system  sensitivity  at  the 

larqer  post-detection  handwidths  needed  to  observe  atmospheric  emission 

fluctuations.  For  most  data  runs,  the  minimum  detectable  temperature 

(  A  T  .  )  of  the  radiometer  was  less  than  I'^K  for  a  6.?F  H?  bandwidth 

rm  n 

^P.lf  sec  intpqration) .  The  actual  AT  .  for  each  data  run  is  noted 

m  n 

on  the  time  hi story/spectral  density  plots  in  Appendix  A. 

Thp  ovprall  nerfomance  of  the  radiometric  measurinq  system  was 
excel Ippt  diirino  the  l^PP  measurements.  System  sensitivity  and 
calibrafion  rpeoatahil ity  were  several  orders  of  magnitude  better  than 
that  obtained  during  thp  107?  and  1970  measurement  programs  at  opL, 
where  serious  radio  froguency  (RFT)  from  the  ambient  electrical 
environment  degraded  system  performance.  >'hile  some  RFI  was  noted 
during  operation  on  thp  Georgia  Tech  campus,  the  effect  upon  the  svstem 
was  not  nearly  as  severe  nor  prolonged  as  that  experienced  at  NRL. 
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R.  FXPrPIMfflTAL  opsFP''ATinM<; 


Ohsfrvat i nns  nF  ;itnospheric  Fluctuations  wore  nadp  undpr  a  varipty 

of  atmospheric  conditions  durinn  the  period  from  February  throuqh 

^’arch  If,  IhBfl.  Appendix  A  lists  the  data  taken  during  this  period  with 

notes  on  weather  conditions  and  the  approximate  AT  .  existinq  For  the 

m  n 

ob'^orvations.  The  nininum  detectable  temperature  was  checked 
periodically  by  observations  on  the  ambient  load,  reference  load  and  hot 
load,  estimating  the  neak-to-peak  noise  N  and  from  this,  calculating  the 

fi 

minimum  detectable  temperature,  aT  .  =  -^  .  (‘^F)  where  SF  is  a  scale 

’  mi  n  6  ' 

Factor  given  in  ^K/inch  of  recorder  paper.  The  weather  conditions, 
listen  in  Appendix  1,  provided  a  large  variety  of  conditions  ranging 

3 

from  clear  skies  te  overcase.  From  low  humidity  (~  0.4  n/m''l  to  high 

g 

humidity  (  ~  IF  g/m")  and  strong  winds.  The  conditions  oF  the  sky 
^clear,  evercast,  ntc.l  can  be  obtained  more  accurately  From  the  notes 
in  Appendix  A,  since  these  data  were  derived  from  observations  made  at 
the  site.  The  data  in  Appendix  B  were  furnished  by  MOAA  at  the  Atlanta 
Airport,  dp  March  ?,  IdBn,  the  weather  conditions  were  very  severe  with 
the  temperature  on  the  order  of  F,  strong  winds  and  snov/  flurries. 

Observations  on  the  ambient  load  were  made  not  only  by  the  load 
switching  scheme  previously  described  [?!]  but  by  placina  a  load 
directly  over  the  antenna  Feed  horn.  This  allowed  an  observation  to  be 
made  to  include  all  losses  within  the  system  (switching,  reflectors, 
horns,  etc.)  except  for  the  antenna  disk.  The  measurements  showed  no 
significant  diFfgrence  from  the  internal  switching  as  far  as  noise  level 
of  the  system  v<as  concerned.  Meither  technique,  switching  or  external 
load  observaFion,  resulted  in  Fluctuation  levels  comparable  to  the  sky 
observations.  These  observations  on  the  ambient  load  indicate  that  RFI 
or  other  external  efFects  did  not  contribute  to  the  large  fluctuations 
observed  during  sky  measurements.  RFI  should  be  the  same  for  the 
external  load  case  and  the  sky  observations. 
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r.  n^TA  ANn  PRFLIMIMARY  analysis 


Thp  ra'ii onrt ri r  Hata  takpn  Hnrim  Lhp  pprioH  from  PS  Pphruarv, 
to  March,  If'RT  arc  nivon  in  Apppnf^icps  T  fAnhient  Load  Pata)  and  P 

(Zenith  Skv  ^ata).  The  hriahtnes<;  tonperature  fluctuations  for  sky 
observation  are  nurh  larper  than  the  fluctuations  obtained  durinp 
observat  i  on<;  op  the  anbient  load.  They  do,  in  fact,  exceed  by  a  larpe 
"’aroin  the  theoretically  expected  values.  Table  1  nives  the  briqhtness 
tefTiperaturo  fluctuations  as  detemined  fron  the  data  of  the  Appendices, 
■'’able  ]A  lists  the  peak-to-neak  tenperature  t)  dct.uati  ons  and  the 
correspondi np  ‘Tpyc.  Spr  nost  cases,  the  anbient  load  observations 
were  on  the  order  of  5  or  less  for  peak-to-peak  fluctuations 

vie'dino  4  1  .P  ”k.  For  the  zenith  sky  neasurenents,  the  peak- 

to-peak  fluctuations  ranqed  between  IS  r  -  30  °K  as  sky  conditions 

varied  considerahl y.  The  resuitinq  ranqed  fron  ?.s  -  S 

^k.  The  rnost  connonly  used  paraneters  durinn  the  observations  were  ?0 
seconds  ^or  the  observation  tine  and  an  inteqration  time  0.16  sec.  In 
most  cases,  the  fluctuation  rate  was  much  faster  than  t  =0.16  sec. 
dn  some  observations,  as  for  Puns  lP-43  on  February  38,  the  temperature 
scales  v/ere  too  low  and  fluctuations  were  excessive  compared  to  the 
laroe  fluctuations  for  zenith  sky  observed  in  other  runs.  It  is  not 
evident  what  the  cause  of  these  chanqes  were,  particularly  since  the 
cbaractori sties  returned  to  typical  values  on  Run  A?  of  that  date. 

The  preliminary  analysis  performed  bv  R.  W.  McMillan  was  intended 
tr>  provide  an  estimate  of  the  fluctuations  which  can  he  expected  for 
reasonable  assumptions  for  atmospheric  parameters.  The  estimate, 
however,  results  in  expected  fluctuations  of  only  a  few  denrees  maximum 
amplitude,  an  order  of  maqnitude  less  than  the  peak-to-peak  amplitudes 
oF  riiirtuation  rhat  have  been  observed.  A  more  accurate  method  of 
determinino  the  meteorolonical  conditions  is  needed  in  order  to  make 
detailed  calculations  of  the  effects. 


Table  lA 

Brightness  Temperature  Fluctuations 


1.)  Ambient  Load  Observations: 
Observation  Time  =  20  sec. 
I  =  0.16  sec 


February 

28,  1980 

Ruji  if 

^^Peak-PeaK 

At 

RMS 

22 

7.5 

1.25 

28 

19 

3.16 

38 

16 

2.6 

57 

6.5 

1.08 

62 

5 

0.83 

67 

6 

1.0 

February 

29,  1980 

74 

8.5 

1.42 

76 

4.3 

0.71 

83 

4.25 

0.708 

84 

5.00 

0.83 

85 

4.5 

0.75 

86 

4.25 

0.708 

110 

9.25 

1.54 

119 

4.4 

0.73 

17 

6.25 

1.04 

18 

6.00 

1.00 

30 


I 

I 

j  March  2,  1980 


34 

6.5 

1.08 

44 

4.5 

0.75 

58 

4.5 

0.75 

Marcli  8,  1980 


19  5.5  0.92 

8  5.0  0.83 


Marc!  i  11,  19 8  0 


8 

5.0 

0.83 

T  =  .5  sec 

13 

5.8 

0.97 

slow  variations 

25 

3.0 

0.50 

T  1.6  sec 

very  slow  variations 

29 

2 

0.33 

42 

8 

1.33 

^  =  0.5  sec 

changed  to  5  second  observation 

52 

21 

3.5 

60 

March  16, 

38 

1980 

6.0 

"  "  "  T  0.016  sec 

78 

4.0 

0.667 

Increased  back  to  20  sec  observation 

99 

3.5 

0.583 

T  =0.16  sec 

115 

2.5 

0.418 

88 

89 

1.0 

0.167 

5  second  observation  low  fluctuations, 
but  rapid  excursions 

2.2 

0.37 

20  sec  observation 

99 

4.1 

0.68 

T  =  0.16  sec 

1 

J 

I 

I 
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Table  IB 

Brightness  Temperature  Fluctuations 


2.)  Zenith  Sky  Observations 

Observation  Time  =  20  sec 
T  =  0.16  sec 


9  24  4 

10  22  3.67 

11  22  3.67 

12  26  4.33 

13  21  3.50 


18-43  see  note  at  bottom  of  Table 


52  23  3.84 

53  35  5.84 

54  30  5.00 

55  30  5.00 

56  30  5.00 

58  31  5.17 

59  28  4.67 

63  26  4.33 

64  22  3.67 

65  23  3.84 

71  27  4.50 

72  26.5  4.41 

73  23.5  3.92 


32 


February  29,  1980 


I 

I 

I 


78 

16 

2.67 

79 

17 

2.84 

80 

17 

2.84 

81 

16 

2.67 

82 

16 

2.67 

87 

16 

2.67 

88 

18 

3.00 

89 

17 

2.84 

90 

18 

3.00 

91 

17 

2.84 

97 

20 

3.33 

98 

18 

3.00 

106 

17 

2.84 

107 

18.5 

3.08 

108 

22 

3.67 

109 

18 

3.00 

116 

16.5 

2.75 

117 

19 

3.17 

118 

17.5 

2.92 

12 

15.5 

2.58 

13 

17 

2.84 

14 

16 

2.67 

15 

17 

2.84 

16 

19.75 

3.29 

! 

I 

I 


I 

I 

I  March  2.  1980 


36 

19.5 

3.25 

37 

20.1 

3.35 

38 

19 

3.17 

39 

20 

3.33 

40 

20 

3.33 

45 

15.5 

2.58 

46 

16.8 

2.80 

47 

15.9 

2.65 

50 

21 

3.50 

51 

19.4 

3.23 

52 

19.25 

3.21 

53 

18.9 

3.15 

54 

20 

3.33 

60 

16.4 

2.73 

61 

17 

2.83 

62 

16.5 

2.75 

63 

15.8 

2.63 

64 

17 

2.83 

65 

21 

3.50 

71 

19 

3.17 

72 

20 

3.33 

73 

18 

3.00 

74 

22.9 

3.82 

75 

17 

2.83 

1 

I 


March  11,  1980 


9 

18 

3 

10 

21 

3.5 

12 

17 

2.83 

18 

9.5 

1.58 

T  =  0.5  sec 

21 

9.5 

1.58 

Change  in  fluctuation  characteristics 

t(  tl  M 

II 

22 

9.0 

1.5 

23 

12 

2.0 

tl  tl  II 

II 

30 

4.0 

0.67 

T  =  1,6  see  " 

II 

31 

5.0 

0.83 

It  tl  II 

It 

32 

7.0 

1.17 

tt  II  tl 

tl 

33 

10.00 

1.67 

It  tl  II 

"  ",  slower 

3A 

5.00 

0.83 

tl  It  tl 

II  II  M 

» 

44 

23 

3.83 

T  =  0.5  sec  " 

"  ",  faster 

45 

23 

3.83 

tl  It  It 

II  II  tl 

» 

46 

24 

4.00 

It  It  tt 

It  It  tl 

> 

47 

24.4 

4.07 

tl  tl  tl 

II  tl  It 

March  16.  1980 

80 

16.7 

2.78 

T  0.16  sec 

81 

17.5 

2.92 

It 

82 

20 

3.33 

M 

83 

20 

3.33 

It 

100 

7.7 

1.28 

60  sec  observation 

with  change  in 

101 

10 

1.67 

characteristics  of 

fluctuation 

102 

8 

1.33 

35 


117 

7.5 

1.25 

T  =  1.6  sec,  60 

sec  observation 

118 

4.0 

0.67 

II  tl 

M  II 

119 

5.5 

0.92 

It  It 

M  M 

82 

10 

1.67 

T  =  5  sec,  500 

sec  observation 

83 

6 

1.0 

tl  (I 

It  II 

84 

4 

.67 

tt  l( 

tl  tl 

101 

16 

2.67 

20  sec  observation,  x  =  0.16  sec 

102 

18 

3.00 

tt  tl 

It  II 

103 

15 

2.50 

It  II 

tl  II 

Note:  For  zenith  observations  18  through  43  on  February  28,  1980,  the 
temperature  scale  was  too  low  and  fluctuations  were  ^  30  K. 
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Uhorpas  the  observations  indicate  that  larqe  fluctuations  occur 
under  the  atmospheric  conditions  that  prevailed  durino  the  February- 
March,  IbSO  period,  a  more  riqorous  analysis  is  necessary.  The 
assumption  that  the  atmosphere  is  stable  durinq  the  observation  period 
must  be  analvzed  further.  V.  E.  Derr  of  NDAA  (Boulder),  in  a  private 
communication  to  R.  W.  McMillan,  has  indicated  that  stronq  fluctuations 
will  be  observed  when  the  atmosphere  is  in  a  state  of  transition,  such 
as  at  the  time  of  formation  of  foqs  and  clouds.  He  has  observed  stronq 
return  usinq  an  R  mm  radar  from  areas  of  clear  sky  near  clouds,  while 
simultaneous  observations  with  a  ruby  lidar  showed  no  returns.  Me 
attributes  these  return  to  refractive  index  inhomoqeneities  caused  by 
water  vaoor.  Such  conditions  could  well  have  existed  durino  the 
measurements  of  this  proqram. 
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1'',  roNri.iisiOM';  AMP  RrmMMFNnA^iOMS 


Observations  performefl  during  the  period  February  -  March  Tf~, 
loon  have  shown  sky  temperature  fluctuations  ranging  from  to 

peak-to-peak  values  with  These  values  correspond 

to  several  sky  conditions  from  clear  to  overcast.  Preliminary 
calculations  indicate  that,  for  clear  weather  conditions,  the  expected 
fluctuations  are  smaller  by  a  Factor  of  approximately  4-7  than  the 
observed  fluctuations.  The  meteorological  conditions  were  such  that 
considerable  instability  could  have  existed  in  the  atmosphere.  The  high 
rate  of  fluctuations  is  also  not  expected.  Some  sky  observations 
appeared  to  have  a  rapid  systematic  variation  characteristic  of 
instability  within  the  radiometer  or  from  external  RFI.  The  lack  of 
such  effects  during  the  ambient  load  observations,  however,  does  not 
support  the  contention  that  the  observed  fluctuation  originated  from 
system/RFl  problems.  Ffforts  to  associate  the  large  fluctuations  with 
sources  other  than  the  atmosphere  have  not  resulted  in  any  conclusions. 
Some  aspects  of  the  radiometer  which  was  used  were  not  desirable  and  do 
present  potential  sources  of  fluctuations.  This  is  particularly  the 
case  for  the  open  structure  beam-waveguide  apparatus  in  the  front-end  of 
the  radiometer.  This  structure  was  employed  for  the  observations  on  the 
10'  dish  at  mpl  but  was  not  needed  for  the  measurements  at  Georgia  Tech. 
It  did  contribute  to  the  degradation  of  the  noise  figure.  For  the  AT^^.^ 
observed  with  the  system,  the  overall  system  noise  figure  was 
approximately  IR  dR.  One  should  expect  at  04  GHz  that  a  NF*:  7-R  dB 
should  be  achievable  with  conventional  radiometric  systems  for  t=  n.lR 
sec.  However,  although  this  structure  did  degrade  the  sensitivity  of 
the  radiometer,  it  cannot  account  for  the  large  variations  (above  the 
minimum  AT)  which  were  observable  when  viewing  the  sky. 
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The  work  of  Kemp  f?!  has  associated  scintillation  with  water  vapor 
fluctuations  in  the  atmosphere.  The  suggestion  of  Herr,  however,  for 
the  transition  case  of  changing  atmospheric  conditions,  would  have  us 
take  into  consideration  refractive  index  changes.  Kemp's  observations 
were  made  during  high  humidity  of  summer  months  whereas  the  measurements 
of  this  program  were  made  for  water  vapor  concentrations  that  varied 
from  n.A  q/m  to  g/m  .  It  is  quite  probable  that  fluctuations 
of  the  magnitude  that  we  have  observed  did  not  originate  from  water 
vapor  fluctuations  alone,  particularly  for  the  case  of  the  low  water 
vapor  concentrations. 

The  rapid  fluctuations  which  have  been  observed  in  these 
experiments  are  also  at  variance  with  theoretical  predictions,  in  that, 
if  the  Taylor  hypothesis  is  assumed  to  apply,  then  the  water  vapor 
distribution  is  assumed  to  drift  through  the  beam  of  the  radiometer  as  a 
result  of  atmospheric  movements  and  the  rate  of  change  in  water  vapor 
distribution  is  considered  to  be  slow  compared  with  the  time  taken  to 
transit  through  the  beam.  The  transit  time  would  have  to  be 
exceptionally  short  for  the  Taylor  hypothesis  to  hold. 

The  data  presented  provides  some  initial  information  on 
fluctuations  at  dHz.  It  is  recommended  that  the  following  tasks  be 
performed : 

1)  Perform  a  more  rigorous  theoretical  study  of  the  fluctuation 
effects  to  define  the  expected  effects  on  a  passive  system  more 
accurately. 

p)  Perform  more  extensive  measurements  under  a  greater  variety  of 
atmospheric  conditions.  It  is  necessary  to  correlate  observations  with 
weather  conditions;  it  is  not  evident  that  all  reports  of  clear  weather 
are  for  conditions  void  of  thin  visually  invisible  clouds. 

'*)  Among  the  observations  to  be  made  are: 
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a.  Simultaneous  radiometric  measurements  at  04  GHz  with 
two  separable  radiometers.  Measurements  as  a  function  of 
the  separation  of  radiometers  are  important  for  determining 
fluctuations  across  the  face  of  a  large  antenna.  Varia¬ 
tion  of  the  separation  of  the  radiometers  should  be  per¬ 
formed,  as  should  interchanging  of  the  position  of  radio¬ 
meters.  The  latter  test  should  remove  the  effect  of  immediate 
surroundings  on  the  apparatus. 

b.  Simultaneous  measurements  at  14n  r,Hz  and  ??n  GHz  to 
examine  fluctuations  as  a  function  of  frequency. 

c.  Measurements  are  needed  for  several  different  weather 
conditions,  actually  for  all  seasons  if  possible. 

d.  On  the  basis  of  Kg-band  radar  observations,  v.R.  Oerr 
of  MOAA  has  recommended  that  simultaneous  radar  and 
radiometry  measurements  be  performed  on  the  same 
propagation  path.  Radar  and  lidar  observations  would 
provide  information  on  particles  not  observable 

visually. 

The  measurements  reported  here  must  be  considered  a  first  effort 
toward  fluctuation  effects  in  the  millimeter  region.  Considerably  more 
data  are  needed.  The  observed  effects  cannot  be  expected  to  occur  under 
all  measurement  conditions  as  the  prevailing  conditions  during  the 
period  that  measurements  were  taken  in  this  work  were  quite  severe  and 
continually  changing  during  observations.  Measurements  as  a  function  of 
seasonal  variations  over  an  extended  period  of  time  are  needed.  A  more 
accurate  method  of  determining  the  conditions  of  clear  sky,  compared  to 
that  of  visual  observation  employed  in  this  program,  is  needed. 

The  work  of  Hill  et  al  ^??1  should  be  investigated  for  extension  to 
this  problem.  values  of  pertinent  parameters  of  their  expressions  are 
needed  for  comparison  with  experiments. 
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APPENDIX  E 


BIBLIOGRAPHY  OF  AFMOSPHERIC  FLUCTUATION  EFFECTS 

A  search  of  papers  relevant  to  this  investigation  has  been 
performed.  The  majority  of  papers  have  been  concerned  with  optical 
propagation  experiments  of  horizontal  one-way  transmission  with  coherent 
sources.  Much  of  the  theoretical  understanding  of  turbulence  effects 
originates  from  these  investigations.  Only  recently  have  experiments 
been  extended  to  millimeter  wavelengths.  Not  many  publications  have 
treated  vertical  radiometric  observations  of  atmospheric  fluctuations, 
and  most  of  these  have  involved  fluctuation  rates  less  than  1  Hz. 
Theory  related  to  vertical  passive  observations  must  also  be  developed 
in  greater  detail.  In  addition  to  turbulence  effects,  cloud  formations 
must  be  included  in  the  formulation  of  sky  brightness  temperatures.  The 
importance  of  wind  effects  is  treated  in  many  publications.  Despite  the 
need  for  advances  in  both  radiometric  observations  and  related  theory, 
most  of  the  papers  of  the  bibliography  present  various  aspects  which  are 
important  to  the  fluctuation  observations. 
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of  Cloud  or  Rains",  IEEE  Trans.  AP,  Vol  AP-25,  No.  5,  pp  650-657,  Sept. 
1977.  Formulation  of  scattering  effects  of  layers  of  clouds  and  rain  on 
down-looking  radiometers.  Derive  radiative  transfer  equations  accounting 
for  polarization  dependence  and  drop  size  distributions.  Models  clouds 
and  solves  resulting  equations  for  brightness  versus  frequency  up  to  300 
GHz. 
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(2)  G.  G.  Haroules,  W.  E.  Brown,  "A  60-GHz  Mul ti -Frequency  Radiometric 
Sensor  for  Detecting  Clear  Air  Turbulence  in  the  Troposphere",  IEEE 
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April  1974.  --  Emission  from  rain  due  to  scattering  effects  at  37,  72  and 
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beam  area  derived. 

(4)  R.  S.  Lawrence,  J.  W.  Strohbehn,  "A  Survey  of  Clear-Air  Propagation 
Effects  Relevant  to  Optical  Communications",  Proc.  IEEE,  Vol.  58,  No. 

10,  pp.  1523-1545,  Oct.  1970.  —  Good  treatment  of  how  inhomogeneous 
regions  are  modeled  in  the  optical  region.  Large  and  small  scale  refractive 
variations  treated. 

(5)  J.  W.  Strohbehn,  "Line-of-Sight  Wave  Propagation  Through  The  Turbulent 
Atmosphere",  Proc.  IEEE,  Vol.  56,  No.  8,  pp.  1301-1317,  August  1968,  -- 
A  review  of  treatments  involving  random  fluctuations  from  optical  to 
millimeter  wavelengths.  Geometrical  models  of  turbulent  mediums  presented. 

(6)  A.  Ishimaru,  "Fluctuations  of  a  Beam  Wave  Propagating  Through  a  Locally 
Homogeneous  Medium",  Radio  Science,  Vol.  4,  No.  4,  pp.  245-305,  April 
1969.  --  Spectral  features  of  index  of  refraction  and  its  effect  on 
fluctuations. 

(7)  H.  J.  Liebe,  J.  0.  Hoppenen,  "Variability  of  EHF  Air  Refractivity  with 
Respect  to  Temperature,  Pressure,  and  Frequency",  IEEE  Trans.  AP,  Vol. 

SP-25,  No.  3,  pp.  336-345,  May  1977  —  Treat  40  to  140  GHz  band  for 
simulated  atmosphere  up  to  40  km.  Present  graphs  of  refractivity  vs. 
frequency,  temperature  and  pressure.  Model  0^  spectral  contributions 
to  refractivity  variations. 
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(8)  C.  B.  Hogge,  R.  R.  Butts,  "Frequency  Spectra  for  the  Geometric 
Representation  of  Wavefront  Distortions  Due  to  Atmospheric  Turbulence", 

IEEE  Trans.  AP-24,  No.  2,  pp.  144-154,  March  1976.  --  Frequency  spectrum 
of  fluctuations  versus  wind  velocity  and  antenna  aperature  derived. 

(9)  L.  Shen,  "Remote  Probing  of  Atmosphere  and  Wind  Velocity  by  Millimeter 

Waves",  IEEE  Trans,  AP,  Vol .  AP-18,  No.  4,  pp.  493-497,  July  1970.  -- 
2 

Effect  of  wind  on  discussed  and  modeled.  Limitations  concerning 
knowledge  of  path  wind  velocities  discussed. 

(10)  MB.  Kanevskii,  "The  Problem  of  the  Influence  of  Absorption  on  Amplitude 
Fluctuations  of  Submillimeter  Radio  Waves  in  the  Atmosphere",  Scientific- 
Research  Radio-Physics  Institute,  Vol.  15,  No.  12,  pp.  1939-1949,  Dec. 

1972.  --  Show  intensity  dependence  of  amplitude  fluctuations  on  wavelength 
due  to  large  and  small  scale  turbulence. 

(11)  A.  0.  Izyumov,  "Amplitude  and  Phase  Fluctuations  of  a  Plane  Monochromatic 
Submillimeter  Wave  in  a  Near-Ground  Layer  of  Moisture-Containing  Turbulent 
Air",  Radio  Eng.  and  Elec.  Phys.,  Vol.  13,  No.  7,  pp.  1009-1013,  1968.  -- 
Increase  in  fluctuations  in  window  regions  compared  to  absorption  line  centers 
discussed. 

(12)  A.  0.  Izyumov,  "Frequency  Spectrum  of  Amplitude  Fluctuations  of  a  Plane 
Electromagnetic  Wave  in  Submillimeter  Range  Propagating  in  a  Surface  Layer 
of  Turbulent  Atmosphere",  Radio  Eng.  and  Elec.  Phys.,  Vol.  14,  No.  10 

pp.  1609-1611,  1969.  --  Presents  frequency  spectrums  for  various  sizes 
of  inhomogeneities  at  300  GHz.  Treatment  has  bearing  on  94  GHz  modeling. 

(13)  L.  A.  Hoffman,  et  al ,  "Propagation  Observations  at  3.2  Millimeters",  Proc. 
IEEE,  Vol.  54,  No.  4,  pp.  449-454,  April  1966.  --  Dry,  wind  atmospheric 
effects  on  signal  scintillation  observed  and  discussed. 
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(14)  T.  Orhaug,  "The  Effect  of  Atmospheric  Radiation  in  the  Microwave 
Region",  Publ ,  of  NRAO,  Vol .  1,  No.  14,  pp.  215-250,  Oct.  1962  - 
Causes  for  fluctuating  component  of  antenna  temperature  discussed. 
Variation  formulated  in  terms  of  fluctuations  in  absorption  coefficient. 

(15)  T.  C.  L.  G.  Sollner,  "Frequency  Spectrum  of  Fluctuation  in  Submillimetre 
Sky  Emission  and  Absorption",  Astron,  Astrophys.,  Vol.  55,  pp.  361-368, 
1977.  -  Dual  beam  astronomy  investigation  of  variations  in  sky  emission. 
Shows  power  spectra  up  to  1  Hz  for  350  g  window  region. 

(16)  N.  D.  Mavrokoukoulakis  et  al ,  "Temporal  Spectra  of  Atmospheric  Amplitude 
Scintillations  at  110  GHz  and  36  GHz",  IEEE  Trans.  AP,  Vol.  AP-26,  No. 

6,  pp.  875-877,  Nov.  1978.  --  Results  of  propagation  experiment  over  a 

4  km  London  path.  Shows  spectral  density  roll-off  at  4.5  Hz  at  110  GHz. 
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Good  theoretical  vs.  experimental  agreement. 

(17)  R.  S.  Cole  et  al,  "The  Effect  of  the  Outer  Scale  of  Turbulence  and 
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Trans.  AP,  Vol  AP-26,  No.  5,  pp.  712-715,  Sept.  1978.  --  Theoretical 
treatment  of  data  from  reference  (16).  Show  how  varying  outer  scale  of 
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no  and  36  GHz",  Elec.  Lett.,  Vol.  13,  No.  7,  pp.  181-182,  March 
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and  temperature  changes  noted. 
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